Mon. Not. R. Astron. Soc. 000,[T]|7]() Printed 5 February 2008 (MN MgX style file v2.2) 



Type II superconductivity and magnetic flux transport in 
neutron stars 

P. B. Jones* 

Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford 0X1 3RH 



ABSTRACT 

The transition to a type II proton superconductor which is believed to occur in a cool- 
ing neutron star is accompanied by changes in the equation of hydrostatic equilibrium 
and by the formation of proton vortices with quantized magnetic flux. Analysis of the 
electron Boltzmann equation for this system and of the proton supercurrent distri- 
bution formed at the transition leads to the derivation of a simple expression for the 
transport velocity of magnetic flux in the liquid interior of a neutron star. This shows 
that flux moves easily as a consequence of the interaction between neutron and proton 
superfluid vortices during intervals of spin-down or spin-up in binary systems. The 
differences between the present analysis and those of previous workers are reviewed 
and an error in the paper of Jones (1991) is corrected. 

Key words: magnetic fields - stars:neutron - pulsars:general 



1 INTRODUCTION 

The origin and modes of evolution of neutron star magnetic 
fields have been topics of interest since the early papers of 
Pacini (1967) and Gold (1968). We refer to Bhattacharya & 
Srinivasan (1995) for a review of the many later papers pub- 
lished on these problems. The difficulty in the interpretation 
of observations on radio pulsars and binary-system neutron 
stars has been that the magnetic flux transport properties 
of both the solid crust and the liquid interior are not well 
known. Empirical deduction of these properties is not fea- 
sible and a priori theoretical input is required. A series of 
papers on the structure of the solid phase in the neutron-drip 
region of the inner crust have shown that it is amorphous, 
heterogeneous in nuclear charge Z, with a high temperature- 
independent resistivity (Jones 2004). But there remain un- 
certainties about the movement of flux in the liquid interior, 
the problem with which this paper is concerned. 

The outer region of the liquid core is believed to con- 
sist of a normal system of ultra-relativistic electrons, a So 
or 3 P? neutron superfluid, and a 1 Sq proton superfluid (see 
Baym & Pethick 1979). Negative /^-mesons, essentially non- 
relativistic, are also present if the electron chemical poten- 
tial exceeds their rest energy. Our assumption, following the 
early paper of Baym, Pethick & Pines (1969), is that the 
protons form a type II superconductor. Type I supercon- 
ductivity will be present in any density region for which 
ft = A/£ < l/\/2, where A is the proton penetration depth 
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and £ the coherence length, but this condition can be satis- 
fied only for very small proton energy gaps. During the early 
stages of neutron-star cooling, the proton phase transition, 
from normal to superconducting Fermi liquid, is accompa- 
nied by the formation of a mixed state. On microscopic 
scales, magnetic flux becomes quantized. The quantum of 
magnetic flux is cj>o = hc/2e = 2.07 x 10 -7 G cm 2 , and it is 
confined to the core of a proton vortex. The highest-density 
regions of the inner crust may consist of low-dimensional 
structures (see Pethick & Ravenhall 1995) in place of spheri- 
cal nuclei. In this case, the boundary of the proton supercon- 
ductor is not well-defined but we shall nevertheless assume 
the existence of a spherical surface within which macroscopic 
supercurrents can flow. The movement of flux across this 
surface depends, of course, on the B(H) characteristic of 
the superconductor and the value of its lower critical field 
H c i (see, for example, Fetter & Hohenberg 1969). 

Apart from the effects of ohmic diffusion and Hall drift 
of the field, we may assume that an approximation to static 
hydrodynamic equilibrium exists at times before the super- 
conducting phase transition. But the transition changes this 
equilibrium because the components of the proton super- 
conductor stress tensor are larger than those of the normal- 
system Maxwell tensor by factors of the order of H c i/B, 
where B is the mean magnetic flux density of the supercon- 
ductor mixed state (Jones 1975, Easson & Pethick 1977). 
The divergence of the superconductor stress tensor is a vol- 
ume force and, with the different transport properties of 
the superconducting system, produces a proton-vortex drift 
velocity. The force has been referred to as a buoyancy force 
and was first estimated by Muslimov & Tsygan (1985). They 
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obtained a drift velocity from its steady-state equilibrium 
with a viscous force derived from the magnetic scattering of 
electrons by an isolated vortex. This specific problem was 
later re-examined by Harvey, Ruderman & Shaham (1986). 
Further work by Jones (1987, 1991) and by Harrison (1991) 
showed that the magnetic viscous force acting on an isolated 
vortex moving relative to the electrons is not relevant to the 
problem of obtaining a drift velocity and that it is necessary 
to consider the interaction of the whole system of vortices 
which is arranged as a two-dimensional lattice. The differ- 
ent ways in which magnetic flux might be expelled, primar- 
ily from normal Fermi systems, were analyzed by Goldre- 
ich & Reisenegger (1992) who emphasized the importance 
of the stable stratification of neutron-star matter and intro- 
duced the concept of ambipolar diffusion, that is, drift of the 
charged components of a plasma relative to the neutral such 
as may occur in the interstellar medium (see, for example, 
Spitzer 1968). The severe constraint imposed on ambipolar 
diffusion by stable stratification was the motivation for the 
work of Ruderman, Zhu & Chen (1998). 

There are several reasons for writing a further paper 
on this problem. The papers cited above, in so far as they 
consider type II superconductivity, are contradictory. There 
is also an error in the paper by Jones (1991). Given the 
quantity of observational data on radio pulsars and on bi- 
nary systems, there is a possibility that a clear solution to 
the problem might reveal interesting information about the 
neutron star interior. The idea that proton vortices move 
easily as a result of interaction with neutron vortices dur- 
ing intervals of spin-down or spin-up is widely assumed in 
studies of X-ray binary systems (see Bhattacharya & Srini- 
vasan 1995) and unambiguous theoretical confirmation of its 
validity, or otherwise, is desirable. 

The approach to the problem made in the present paper 
uses the fact that, at the temperature considered here, all su- 
perfluid quasiparticles except those localized in vortex cores 
have negligibly small number densities. Thus it is possible to 
write down a Boltzmann equation for the electrons (muons) 
interacting only with muons (electrons) and with the pro- 
ton vortices and obtain a steady-state solution (Section 2). 
The screening of the lepton current in the superconductor 
and the proton force-balance give further equations connect- 
ing drift velocities and chemical potential gradients for lep- 
tons and protons (Section 3). A simple and unambiguous 
result is obtained for the drift velocity of the proton vor- 
tices (equation 22) and hence for the expulsion of flux from 
the interior of the neutron star. It is the principal result of 
this paper, and is shown to be unaffected by possible lep- 
ton interactions with other degrees of freedom. In Section 
4, we attempt to analyze the differences between this paper 
and the work of each of the sets of authors cited above. The 
scope of the paper is purely technical and it does not con- 
sider observational evidence relevant to magnetic flux evolu- 
tion. It excludes the possibility of type I superconductivity 
induced by interaction between neutron and proton super- 
fluids (Buckley, Metlitski & Zhitnitsky 2004), also systems 
of several distinct superconductors which may be present at 
very high matter densities in the inner liquid core. Unless 
otherwise stated, all quantities are defined in a coordinate 
system corotating with the solid crust of the star and having 
angular velocity ft. 



2 THE LEPTON BOLTZMANN EQUATION 

The reference system we shall use within the chosen coor- 
dinates is the steady state with zero vortex drift velocity, 
vl = 0. A local definition of electron chemical potential 
is used, ^ = tfe, where ep e is the Fermi energy, includ- 
ing rest energy. For present purposes, this is more conve- 
nient than the global definition, (i e — (e_Fe — e<f>) y^goo, in 
terms of the time-like component of the metric tensor and 
the gravitationally-induced electric potential <f>. This lat- 
ter definition was favoured by Harrison (1991); its equi- 
librium value m e c 2 y/goo s , defined at the stellar surface, is 
constant throughout the star. In the liquid core, eFe can 
exceed the muon rest energy so that a muon chemical po- 
tential is also defined. At these energies, lepton transport 
can be considered under the assumption of non-quantizing 
fields (see Potekhin 1999). The classical leptonic orbits in 
a plane perpendicular to the proton vortices are irregular 
polygons whose vertices represent scattering by the micro- 
scopic magnetic flux density B localized within the vortex 
cores. To be specific, we consider first the electron compo- 
nent. Its orbit size can be specified by the orbit radius in 
the spatially-averaged magnetic flux density B = (B) . This 
is r B = £Fe/eB = 3.3 x lO -7 ^ 1 cm, for e Fe = 100 MeV, 
where B12 is the magnetic flux density in units of 10 12 G. 
It is several orders of magnitude larger than the intervortex 
spacing in the triangular lattice, d — 4.9 x 10~ 10 -B 12 cm. 
For such a system, it is possible to define, within a small el- 
ement of phase space, a spatially-dependent electron Fermi 
distribution function which satisfies a Boltzmann equation 
(see, for example, Pines & Nozieres 1966). In the steady 
state, with vl = and isotropic distribution function n\, 
this can be expressed as 

Vk • Vn" - e (s + iv k x B^ • V k n2 = 

-T° k + e (ivk x (B - B)) • V k n£, (1) 

with electron velocity v k for momentum k. The conservative 
field £ cancels the Vn° term in the equation. Motion of the 
vortex lattice produces an induction field E = — (1/c)vixB 
and changes the electron distribution function to 

nk = n° k + /3n°(l - n°)<5ju e + 5n k , (2) 

in which the first incremental term is isotropic, with /3 _1 = 
fcsT, and represents the effect of a chemical potential change 
5 fi e = SeFe — e5(f>. The second term is anisotropic. The mod- 
ified Boltzmann equation is then, 

v k ■ Vri k - e (e + E + iv k x • V k n k = 

-r k + e (e- E+ iv k x (B -B)) • V k n k (3) 

Both equations have been expressed in terms of the 
spatially-averaged fields. The reason for this is that the dif- 
ference between the irregular polygons and the circular or- 
bits of the spatially-averaged field can be thought of as the 
consequence of a scattering process. Therefore, the difference 
terms producing it are placed on the right-hand side of each 
equation with the collision integrals F k and F k . These are 
derived from processes other than magnetic scattering and 
include electromagnetic scattering by muons (if present) and 
by vortex-core quasiparticles. But these Boltzmann equa- 
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tions are for a distribution function n k (r) defined within a 
small element of phase space. Spatial integration over that 
small element leaves the terms in them unchanged with the 
exception of the difference terms, whose integrals clearly 
converge to zero as the linear dimension of the element be- 
comes large compared with the intervortex spacing d. This 
serves to establish the intuitively obvious result (Harrison 
1991, Jones 1991) that, given the condition rg >• d, the mi- 
croscopic fields E and B can be replaced without error by 
the spatial averages E and B. 

With neglect of some terms of the second order of small- 
ness, the incremental change in distribution function (5n k 
satisfies 



- eSE • V k n£ (vk x B) • V k <5n k = T k - T k , 

c 

where, 



(4) 



5£ = E+ -V5/tt e . (5) 

In order to write down an explicit form for the collision inte- 
gral difference, it is necessary to specify each set of degrees 
of freedom with which the electrons interact. At tempera- 
tures well below the superconducting transition, superfluid 
quasiparticle densities are negligible except for those local- 
ized in vortex cores. The rest frame for these excitations is 
therefore that of the vortices, with velocity v L in the ref- 
erence frame corotating with the solid crust in which our 
quantities are defined. Thus interaction with them causes 
<5n k to relax to isotropy in the rest frame of the vortex lat- 
tice with relaxation time r*. There is also interaction with 
muons (if present) and possibly with other excitations, such 
as zero sound phonons or other collective modes, or with su- 
perfluid continuum quasiparticles whose density will be sig- 
nificant near the transition temperature. In order to include 
them, we assume for the muons a natural rest frame with 
velocity v M , in which their thermal-equilibrium distribution 
function is isotropic, and a separate relaxation time r^. For 
the unspecified excitations, we assume a rest-frame velocity 
v„ and relaxation time r^, with both quantities treated as 
unknowns. Thus the collision integral difference is 

r k -T k = — \ (<5n k - /3n°(l - n°)k • vl) 
-\ (Sn k - (3n° k (l - n°)k • v M ) 

-\ (<5n k - /3n°(l - n°)k- v u ) ■ (6) 

Equations (4) and (6) are satisfied by, 
dn° 

<5n k = er„ e — -^v k ■ A, (7) 
oe k 

where A is related to the electron current density, 

J e = -N e ev e = — % / d 3 fcv k n k = ol A, (8) 
( i7I T J 

and N e is the electron number density. Equations (8) define 
the electron drift velocity v e . The resistivities are 
defined by 



(9) 



electrons interact. The relation between the drift velocities 
Ve l( j, u and v L derived from equations (4) and (6) can be 
expressed in the form, 

1 _ , N e e N e e 

-V8fJ, e ~ — (V M - V e ) - — (V U - V e ) = 

^ (VL - V e ) + -(Vi- V e ) X B. (10) 
(J v C 

This unremarkable result is no more than a force-balance 
equation relating the force on the electrons from the vortex 
lattice (right-hand side) with the forces exerted by the elec- 
trons on the chemical potential gradient, on the muons, and 
on the unspecified system of excitations whose rest-frame 
has velocity v„. A similar equation exists for the muons, 

-V<5/i M — (V e - V 

6 (7 e 



V U 



(VL - V„) + - (VL - V M ) X B. 



(11) 



The rate at which the lattice, the muons and the unspecified 
excitations do work on unit volume of the electrons less the 
rate at which the electron system does work in moving on 
the chemical potential gradient is, 



N, 



( N ee . 1 . , \ 

eV L ■ I — r (VL - Ve) + - (VL - Ve) X B j 
-N e V e ■ V(5/i e + V V (V„ - Ve) 



, (N e e) 2 

~\ (V„ - V e ) • V u 

— — (Vi-Ve) +— — ( 



(12) 



which is the rate of energy dissipation per unit volume at- 
tributable to v e . A similar relation exists for the muons. 

The analogous force-balance equation for the proton su- 
perfluid continuum can be written down immediately be- 
cause the continuum quasiparticle number density is negli- 
gibly small. It relates the chemical potential gradient with 
the Magnus force per unit volume, 



N p \7Sfip = (v p o — vl) x B, 



(13) 



where the subscript denotes the system with which the 



in which v p o is the proton drift velocity. Equations (10), (11) 
and (13) are the principal results of this Section but contain 
a total of 7 variables excluding v„. The further equations 
required can be obtained by examining the effect of the nor- 
mal to superconducting transition on the various current 
densities in the crust and liquid interior. 



3 STEADY-STATE EQUILIBRIUM OF THE 
MAGNETIC FLUX DISTRIBUTION 

Before the transition to superconductivity, the magnetic flux 
and lepton current density are slowly-varying functions of 
position satisfying Ampere's theorem. If the transition is 
to a type II superconductor, it is accompanied on a micro- 
scopic scale by the formation of proton vortices. These have 
quantized magnetic flux <p , locally collinear with their axes, 
supported by a circulating supercurrent distribution whose 
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density decreases exponentially with radius on the scale of 
the penetration depth A. If, on macroscopic scales, the vor- 
tex number density N v is such that the spatially-averaged 
magnetic flux density is unchanged, so that N v <p — B, the 
spatial average of the individual vortex supercurrent distri- 
butions J pa must satisfy the relation, 



= j e +j p 



(14) 



The pre-existing lepton current density must be screened 
out by a supercurrent density J p0 which, in the body of the 
superconductor, has significant variation only over lengths 
many orders of magnitude larger than A, 



j p0 + r + r 



o. 



(15) 



This ensures that Ampere's theorem is satisfied at all points 
in space (see Jones 1991) and, as a consequence of super- 
conductivity, is analogous with the formation of surface cur- 
rent sheets in the Meissner effect. An additional condition 
expresses the fact that protons do not cross the supercon- 
ductor boundary, 



(16) 



The effect of equation (16) is that the supercurrent distribu- 
tion J p0 defined by equation (15) must be associated with a 
return current sheet at the boundary separating the super- 
conductor from the normal solid. This current sheet, which 
has negligible kinetic energy, also maintains the mixed state 
of the superconductor by excluding the continuous magnetic 
flux density of the normal system. 



The neutrality condition N e + 7V M 



N„ can be re- 



garded as exact and equation (15) gives the condition 
N e v e + iV M v M = NpVpo which is satisfied at all points in 
the superconductor by the lepton and proton drift veloc- 
ities. The final equation is for steady-state equilibrium of 
the vortex lattice under what is referred to as the buoyancy 
force. The argument for the existence of this quantity is as 
follows. The normal-state hydrostatic equilibrium before the 
transition is given by the equation, 

w ^£" + i ((je+jM)xB) =°' (i7) 

in terms of the matter density p, the gravitational accelera- 
tion g, and the zero-field pressure P°. The third term is the 
divergence of the Maxwell stress tensor. In hydrostatic equi- 
librium, its curl is almost exactly a non-radial vector. (This 
condition does not apply in the crust owing to the pres- 
ence there of a further term, derived from the solid stress 
tensor, which allows more general components in the diver- 
gence of the Maxwell tensor and the resultant phenomenon 
of Hall drift. In the liquid, V x ((J e + J M ) x B) can pro- 
duce Hall drift only of non-radial components of B. Here, for 
brevity, components in the electromagnetic current density 
other than J e ' M have been neglected.) At the transition, the 
magnetic part of the superconductor stress tensor replaces 
the Maxwell tensor. Following Easson & Pethick (1977), we 
express this in the form, 



— —P s 8ij + -^—HiBj, 



^3 



(18) 



which is a symmetric tensor because B and H arc locally 
parallel vectors. Both isotropic and anisotropic components 
are larger than the corresponding Maxwell components by 
factors of the order of H c i/B 3> 1. The new hydrostatic 
equilibrium is given by 

dP° 



P9i ~ 



dxi 



dTf, 
dxj 



(19) 



where fy is the force per unit volume arising from inter- 
action between neutron and proton vortices (Sauls 1989), 
an additional effect which is not included in equation (18). 
(Strictly, the existence of fv =fc depends on rotation of 
the neutron superfluid with angular velocity Cl n ^ SI.) The 
derivative of the isotropic component of the stress tensor 
has been referred to as a buoyancy force, but the implicit 
neglect of the anisotropic component means that this de- 
scription is not necessarily apt. Easson & Pethick obtained 
an expression (equation 17 of their paper) for the isotropic 
component of the stress tensor valid for any n > l/v2 in 
the limit B <C H c i. But they note that it can be estimated 
reliably only in the extreme type II limit in which k 3> 1. 
Thus for more general values of k, it is not obvious that the 
term leads to a buoyancy force, rather than the reverse. In 
the limit B <C H c \, the magnetic field H w H c \ and we can 
assume that the spatial derivatives of its magnitude (though 
not direction) are much smaller than those of B. With ne- 
glect of these terms, and in the further limit k >• 1, the 
divergence of the stress tensor given by equation (18) is, 



9T% 
dxj 



1 ((V x B) x Hdjj + -!- (V • H c i) B i: 

47T 



47T 



(20) 



in which H c i has magnitude H c \ and is everywhere parallel 
with B. This replaces the third term in equation (19). The 
sum of the third and fourth terms in equation (19) is the 
total magnetic flux-dependent force per unit volume (defined 
as fs) and acts on the vortex lattice. Its direction is not 
necessarily the radial direction expected of a buoyancy force 
and it is analogous with the third term in equation (17). 

The lattice drift velocity vz, is then defined by equat- 
ing fs with the sum of the force on the electron and muon 
systems and the Magnus force on the superfluid continuum, 

„ {N e ef N e e . 

ffl - ^~ (V L - V e ) (VL - V e ) X B 

<J v c 



VV) - (VL - , 

N P e i ^ -d 

— ^— (v p o -vl)xB 



x B 



0, 



which gives an extremely simple final result, 



_ (TffiVgVe + agiV>M 



(N e + N^e- 



where, 



(iV e + iV M ) 2 a>£ 



(21) 



(22) 



(23) 



It is interesting that this result is quite independent of the 
existence of the resistivities derived from lepton interaction 
with other degrees of freedom whose natural rest frame is 
not that of the vortex lattice. Equations (10), (11), (13), 
(15) and (21) give vl and a complete description of the sys- 
tem in terms of the independent variables v e , v M and fs, 
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which all depend on the properties of the normal-state B- 
distribution that existed before the superconducting transi- 
tion. Scattering transition rates for e — p and /i — p at the 
Fermi surfaces, though of the same order of magnitude, are 
not identical so that we anticipate v e 7^ v M and retain both 
as independent variables in equation (22). Of the remaining 
variables v p o and V5/tt e , M , P , the increments in chemical po- 
tential, 8fj. e (SN e ,S<f>), Sfi^SN^jSip) and 5n P (SN p , —S4>), are 
functions of just 3 further variables SN e , <5iV M and 5<fi. The 
values of a that are immediately relevant are those at tem- 
peratures below, but within an order of magnitude of the 
superconducting transition temperature. We refer to Jones 
(1991) for the order of magnitude of crj and for a brief discus- 
sion of the factors determining its temperature dependence. 
But it is possible to assert that as the star cools, it increases 
at least as rapidly as T~ 2 , where T is the temperature. Thus 
at a typical transition temperature of 3 x 10 9 K, the order of 
magnitude is a ~ lO 29 !?^ 1 s _1 . A force with order of mag- 
nitude /s w 10 20 Bi2 dyne cm -3 then leads, from equation 
(22), to a velocity vl ~ 4 x 1CF 7 cm s _1 . It is evident that a, 
and therefore vl both increase rapidly as the star cools. Ex- 
amination of the orders of magnitude present in equations 
(5), (10) and (11) shows that the chemical potential gradi- 
ents can become quite large. But they are almost exactly 
cancelled by the induction field. Our conclusion is that, for 
any plausible value of /s, the post-transition movement of 
magnetic flux is fast in comparison with, for example, radio 
pulsar spin-down or binary-system evolutionary time-scales. 



4 CONCLUSIONS 

Before proceeding to a comparison of the results obtained 
here with those of earlier papers, it is worth considering 
two questions about flux movement. Konenkov & Geppert 
(2001) observed that proton vortices, given a buoyancy force, 
always move toward the crust by sliding (with some dissipa- 
tion) along the rectilinear neutron vortices of the rotating 
neutron superfluid. (The only exceptional case is that in 
which proton vortices, within a small element of solid an- 
gle, are approximately parallel with the neutron vortices.) 
This interesting possibility clearly depends on fs having a 
component, parallel with the neutron vortices and the spin 
vector ft n of the neutron superfluid, with the appropriate 
sign. But we have seen, both from the discussion given by 
Easson & Pethick and by examination of equation (20), that 
this sign must depend on the form of the flux distribution 
prior to the superconducting transition. For example, the 
force terms in equation (20) vanish in the case of a uniform 
B, giving fs = fv, which is almost exactly perpendicular 
to tt n . The second question concerns whether or not proton 
vortices can be moved inward or outward, relative to the 
rotation axis of the star, by interaction with neutron vor- 
tices during intervals of spin-up or spin-down. The neutron 
vortices at a distance r± from the rotation axis move with a 
radial velocity component v± — —r±0,„/20, n , where Sl„ is 
the neutron superfluid angular velocity at r± . From equation 
(22), vl can have a component of this size provided, 



where m p is the proton mass and fv is the maximum force 
that unit length of neutron vortex can exert on a lattice 
of proton vortices without intersection. (Here, we assume 
Sl n = SI, the angular velocity of the star.) It is easy to 
confirm that this is clearly satisfied in the case of the ob- 
served radio pulsars by evaluation for the Crab, which gives 
/yer > 8 x 10 44 dyne cm" 1 s" 1 . A value fv ~ 10 14 dyne 
cm -1 is a conservative assumption (see Jones 1991) which, 
with the estimated a > lO 32 ^^ 1 s _1 for iV M <JC N e , shows 
that the condition is satisfied. However, much higher spin- 
down rates exist during the propeller phase of binary sys- 
tems. The model of Urpin, Geppert & Konenkov (1997) gives 



fl = -(GMf^R^B^M^r 1 , 



(25) 



where M, R and I are, respectively, the neutron star mass, 
radius and moment of inertia. For typical values (M — 
1.4M Q , R = 1.2 x 10 e cm, I = 10 45 g cm 2 and an accretion 



rate on to the Alfven surface of M 
inequality (24) gives, 

fva > 3 x 10 46 B 2 / 7 n -2 , 
that is, fv<7 > 8 x 10 48 dyne cm" 



1O" 1O M yr" 1 ) the 



(26) 

for the rotation 

periods of 10 2 s that are observed at the end of the spin- 
down phase of binary evolution (see the review of Verbunt 
& van den Heuvel 1995). Satisfaction of this condition is 
much more problematic. 

Previous papers concerned with the flux transport ve- 
locity have reached differing conclusions. In their seminal 
paper on type II superconductivity in neutron stars, Baym, 
Pethick & Pines assumed that the viscous force on unit 
length of moving proton vortex would arise from a dissi- 
pation rate of the order of 



t>oH c 2<7 



vl 



(27) 



analogous with the theory of Bardeen & Stephen (1965). In 
equilibrium with the force per unit length fB(f>o/B derived 
from equation (20), it would give extremely small velocities 



(vl 



10" 



cm s ). In this expression, E is the micro- 



fvo > 



nr± I Q 1 he 2 (N e + iV M ) 2 
4m„fi 2 



(24) 



scopic induction field, a is the normal-state conductivity, 
and H C 2 the upper critical field of the superconductor. It 
assumes, in a laboratory superconductor, relaxation of the 
vortex-core electron distribution function to isotropy in the 
frame of the ion lattice. But no equivalent interaction ex- 
ists in a proton superconductor and the theory is therefore 
inapplicable (Jones 1987, 1991: Harrison 1991). 

Muslimov & Tsygan (1985) considered vl to be deter- 
mined by the steady-state equilibrium between a buoyancy 
force (the derivative of the isotropic part of equation 18) and 
a viscous force derived from magnetic scattering of electrons 
by an isolated moving vortex (see also Harvey, Ruderman & 
Shaham 1986). These calculations gave v L ~ 10~ 10 - 10~ 8 
cm s _1 but were subject to the objection that the viscous 
force used was not correct for the problem concerned, which 
was the motion of the whole triangular lattice, with fairly 
small intervortex spacing d — 4.9 x 10 -10 B^^ 2 cm, through 
the electron gas (Jones 1987, 1991; Harrison 1991). These 
authors considered the coupling of the lattice with the elec- 
trons and recognized that, given the limited interactions oc- 
curring in a superfluid system well below the critical tem- 
perature, the electron distribution function must relax to 
isotropy in the rest frame of the lattice. Jones (1987) ne- 
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glected the supercurrent screening condition (equation 15) 
and is therefore seriously in error, as is Jones (1991) who 
obtained a result equivalent to equation (10) but then incor- 
rectly identified the quantity J e -5£ as the rate of dissipation. 
The final conclusions of both these papers must therefore be 
disregarded. Harrison (1991) considered the gravity-induced 
electric field (the term £o of equation 3) which must be 
present in both normal and superconducting systems. His 
treatment of the electron-lattice interaction is based on as- 
sumptions similar to those of Jones (1991). It includes a 
Lorentz force on the lattice, but differs from Jones (1987, 
1991) and the present paper in not considering the proton 
superfluid Magnus force produced by the difference v p o — vl • 
This is an effect arising from classical hydrodynamics (see 
Nozieres and Vinen 1966, Jones 1991) and there would ap- 
pear to be no doubt of its existence. Harrison also observed 
that the buoyancy force must be included in the equation 
for hydrodynamic equilibrium, as in equation (19). His final 
conclusion is that vl <§C 10~ 16 cm s _1 . 

The paper by Goldreich & Reisenegger (1992) is con- 
cerned primarily with matter whose composition is defined 
as being chemically homogeneous (neutrons and charged 
components limited to electrons and protons), without su- 
perfluidity. It introduces to neutron star physics the concept 
of ambipolar diffusion (see also Haensel, Urpin & Yakovlev 
1990). This is movement of the charged components and 
magnetic flux relative to the neutral part of the system. It 
would also be of relevance to the intermediate state in proton 
type I superconductivity if that were present in the neutron 
star interior. These authors observe that, in analyzing such 
motion, it is essential to divide particle flux vectors NiVi 
into solenoidal and irrotational components. Then solenoidal 
motion, for chemically homogeneous systems, is limited only 
by a viscous force. Its velocity is v s ambip oc r pn , where r pn 
is derived from p — n scattering and is, in the normal sys- 
tem, the principal collision relaxation time. (The distinction 
is that solenoidal components do not change local number 
densities, whereas irrotational components produce chemical 
potential gradients and an imbalance Sp = Sp p + <5/i e — 5p n , 
which is the basis for the relationship between stability and 
stratification emphasized by these authors.) 

The structure of the liquid core may be more complex 
than we have assumed in the present paper. Even if type II 
superconductivity were present in the outer region, protons 
in the inner core might be of type I. Type II flux quanta <f>o, 
on entering the inner region, would merge to form an inter- 
mediate state of the type I superconductor in which mag- 
netic flux is confined to filaments of normal proton Fermi 
liquid with macroscopic cross-sectional area. The relaxation 
time Tp„ derived for proton interaction with superfluid neu- 
trons (well below the critical temperature) is extremely long 
so it might be thought that there would be rapid ambipo- 
lar diffusion of these filaments under stresses generated by 
flux movement in the outer (type II) region. But this is not 
necessarily so owing to the presence, in many equations of 
state, of muons at inner-region matter densities. The mat- 
ter is then chemically inhomogeneous in the sense considered 
by Goldreich & Reisenegger. Radial motion would produce 
a large chemical potential difference 5p e — Sp^ because the 
muons are non-relativistic. The consequence, following Gol- 
dreich & Reisenegger, is that even solenoidal particle fluxes 
are strongly inhibited by the stratification and stability con- 



dition so that flux movement out of the type I region would 
be slow, depending on weak-interaction transitions to re- 
move the imbalance. The rates for these are strongly sup- 
pressed by the proton (and possibly also the neutron) energy 
gap; also by the requirement for energy-momentum conser- 
vation in direct p ^ e transitions on the Fermi surfaces with 
neutrino-pair creation. 

We emphasize that the flux velocity vl given by equa- 
tion (22) can be of an order of magnitude different from the 
individual drift velocities v e ,^, p o- In this respect, the motion 
differs from the ambipolar diffusion described by Goldre- 
ich & Reisenegger. On the other hand, there is a similar- 
ity between equation (22) and ambipolar diffusion which is 
worth noting in relation to the proton-vortex drift veloc- 
ity obtained by Ruderman, Zhu & Chen (1998; equations 
10 and 14 of that paper) which is inversely proportional to 
an effective conductivity a. It is defined by dissipation in 
an induction field rather than by a viscous force. The ex- 
planation for this appears to be that, in considering proton 
vortex motion, these authors have not included both the 
Magnus force and the screening condition given by equa- 
tion (15) of the present paper. This appears to be the origin 
of the disagreement with our drift velocity wl- Also, the 
effective conductivity of Ruderman, Zhu & Chen is deter- 
mined principally by the component of the magnetic force 
on an isolated moving proton vortex that is antiparallel with 
its velocity, and for that reason is referred to as the mag- 
netic viscous or drag force. The view of the present paper 
(see also Jones 1991) is that this scattering is the process 
that produces the polygonal electron orbits and, of course, 
is the source of their irregularity. The discussion which im- 
mediately follows equation (3) expresses our view that this 
process, represented in equations (1) and (3) by the terms 
containing B — B, does not contribute to the Boltzmann 
collision integral. 

The work of this paper has been limited to density re- 
gions containing only electrons, muons, protons and neu- 
trons. Most equations of state have hyperon thresholds that 
are likely to be exceeded in a typical 1.4M Q neutron star. 
Thus a complete discussion of flux transport would need to 
consider not merely the possibility of type I proton super- 
conductivity in the higher- density regions, but the physics 
of systems with two or more superconducting baryonic com- 
ponents. 
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